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Abstract 

A generalization of the Kontsevich Airy-model allows one to compute the 
intersection numbers of the moduli space of p-spin curves. These models are 
deduced from averages of characteristic polynomials over Gaussian ensembles 
of random matrices in an external matrix source. After use of a duality, and 
of an appropriate tuning of the source, we obtain in a double scaling limit 
these intersection numbers as polynomials in p. One can then take the limit 
p — — 1 which yields a matrix model for orbifold Euler characteristics. The 
generalization to a time-dependent matrix model, which is equivalent to a 
two-matrix model, may be treated along the same lines ; it also yields a 
logarithmic potential with additional vertices for general p. 
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1 Introduction 



Many topological invariants have been computed from matrix models of mod- 
uli spaces. The well-known Kontscvich's Airy matrix model [1] gives the in- 
tersection numbers of the moduli spaces of curves, which was also studied by 
the double scaling limit of one matrix model [2] . The Euler characteristics 
of orbifolds have been computed from the Penner model [3] . The intersection 
numbers for the p-spin curves are obtained from the generalized Kontsevich 
model [4]. These matrix models all provide explicit results for the intersection 
numbers. 

In this article, we discuss the above models, which lead to three kinds of 
matrix models, in a unified way. Our formulation starts from simple Gaussian 
matrix models with an external matrix source. In recent articles we have 
already considered the average characteristic polynomials in these Gaussian 
ensembles, and derived, through a duality relation and the replica method, 
the intersection numbers of p-spin curves [5, 6, 7]. This duality relates the 
average of the product of k characteristic polynomials ior N x N random 
matrices M, to the average of the product of N characteristic polynomials 
over k X k Gaussian random matrices B. In the large N limit, the matrix 
model for B reduces to the higher Airy matrix [6, 7] for the intersection 
numbers of spin curves studied by Witten[4, 8]. 

This duahty allows one to compute the intersection numbers for the spin 
moduli spaces with n-marked points and genus g, from an n-point correlation 
function f/(si, s„) of Gaussian random matrices in a scaling limit near 
critical edges [9, 10]. The basic steps are recalled in section 2. 

In this article we first compute explicitly the intersection numbers of 
moduli space of p-spin curves with one marked point, for arbitrary values of 
p, as polynomials in p. We have obtained earlier the intersection numbers 
for p—2,3 and 4 explicitly, but we discuss here arbitrary p. This allows us 
to consider continuations in p ; in particular the limit p — — 1 exhibits 
an interesting relation between the intersection numbers, (r-class) and the 
orbifold Euler characteristics x(Mg_i) = ^(1 — 2g) (C(a^) is Riemann zeta 
function) [3, 11]. In section 3 we derive these numbers for surfaces with one 
marked point. 

In section 4 we show that the intersection numbers with n-markcd points 
for p-spin curves are also obtained easily from the integral representation 
of U{si, ...,Sn) at the critical values tuned through an appropriate external 
source. We evaluate the case of two marked points for genus one (g=l) and 
any p. The generating function U is given for three and four marked points 
in Appendices A and B. Our results through this generating function U are 
consistent with the previous recursive results of Virasoro equations. We find 
that the ring structure of the primary fields, for genus zero, is deduced from 
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these n-point correlation functions U for arbitrary p. This shows that the 
random matrix theory with external source near critical edges, has a structure 
of a minimal N — 2 superconformal field theory with Lie algebra Ap_i . 

It has been conjectured by Wittcn that the free energy F which gener- 
ates the intersection numbers of the moduli space of p-spin curves satisfies 
a Gelfand-Dikii hierarchy [4, 12]. Wc show here that the intersection num- 
bers, computed from the integral representation of U{si, Sn), do satisfy 
Gelfand-Dikii equations. We present in Appendix C, this Gelfand-Dikii hier- 
archy equations and the construction of the super potential for the primary 
fields. We also note that, with respect to Witten's conjecture, that the defini- 
tion of intersection numbers as an integral over the compactified moduli space 
Mg^n, is similar in structure to the integral representation of U{si, 

In section 5, begins a second part of the paper devoted to the the time 
dependent Gaussian matrix model for which we extend our previous work 
on duality [5, 6, 7]. The time-dependent model, a matrix quantum mechan- 
ics of harmonic oscillators, reduces easily, for a Gaussian distribution, to an 
equivalent two-matrix model. Again one may derive (section 6) a dual model 
in the presence of a matrix source. We then obtain, with an appropriate 
tuning of the source, the two-matrix equivalents of the Kontsevich plus Fan- 
ner models for matter with central charge c = 1 for the p = 2 case. For 
p > 2, additional terms are present with respect to the c=l matrix model for 
tachyon correlators [14, 15]. 



2 Replica and duality for the one-matrix model 

Let us first summarize the steps which one uses to obtain Airy and higher 
Airy matrix models from the Gaussian one- matrix model in a source, followed 
by the computation of the intersection numbers through the replica method 
[5, 6, 7]. 

The m-point correlation functions of the eigenvalues in the Gaussian 
unitary ensemble are conveniently deduced from their Fourier transforms 
U{si, Sm), defined as 

U{si, S2, ...,Sm) = < tre'^i^^tre*-^^^ ■ ■ ■ tre"'"^^ > 

/m m 
Yl dXie^'''^' < l[trS{Xj -M)> (1) 
1=1 1 

where si = iti ; M is an N x N Hermit ian random matrix. The bracket 
stands for averages with the Gaussian probability measure 

< X J dMe->^'+''''^^X{M), (2) 
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A is an X external Hermitian source matrix. We may assume that this 
matrix A is diagonal with eigenvalues aj. We consider later the external 
source A with (p — 1) distinct eigenvalues, each of them being times 
degenerate. This parameter p is crucial in this paper. 

Let us consider first the one point case (m=l), namely U{s). The replica 
limit k ^ for < tr6{X — M) > relies on the identity 

lhnl([det(A-I-Mr-l) 

= trlog(A • I - M) (3) 

Taking a derivative with respect to A, it yields the density of states tr5(A — 
M) as the imaginary part of the resolvent when the imaginary part of A 
goes to zero. Thus U{s) is expressed in terms of products of characteristic 
polynomials, 

U{s) ^}^Jrj d\e'^§y^ < li det(A« - M) >a U.=a (4) 

a=l 

We have introduced a replica symmetry breaking by taking k distinct A^ 
{a — 1, k) [16], in order to use the duality formula [6, 17], 

I k 1 ^ 

— < n det(Aa - M) >A^ ^ < n det(% - iB) >a . (5) 

^0 a=l ^0 j=l 

where A = diag{Xi, A„), B is an k x k Hermitian random matrix, and Zq 
and Zq are normalization constants. The probability distribution for B in 
the right hand side is 

The simplest case consists of taking an external source multiple of the iden- 
tity, namely aj — 1, j = 1, ■ ■ ■ N. The effect of this source is simply to shift 
all the eigenvalues of M by one : the left edge of Wigner's semi-circle is now 
at the origin. In the large N limit, after exponentiation, expansion of the 
integrand in powers of B, cancellation of the trB^ terms, we obtain 

Z = lim < [det(l - iB)f >A= / ^5g^trB3+iiVtrB(A-l) 
N—*oo J 

which is Kontsevich Airy matrix-model. (We explore here a "double scahng 
limit", namely the vicinity of the origin, in which the (A^ — 1) are of order 

regime NtvB'' is negligible for / > 4). This 
Airy matrix model has an expansion in terms of the moduli parameters tm, 

1 1 
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where C is a normalization constant, to be determined later. When we set 
the k distinct Aq to a common value A, we have tm = C ^2m+i ■ Then, in the 
limit A; — > 0, only single r's appear, since the parameter tm is proportional 
to k. Thus the zcro-rcplica limit A; — > yields the intersection numbers with 
one-marked point (one tm) : 

U{s) = j d\e'^—[l+ < n > ti+ < r4 > t4 + • • •] (10) 

where < ri >= ^, < T4 >= '(2^p^_- As found in [5], U{s) is obtained, after 
approriate normalization, as 

U{s) = ^^ei (11) 

which gives the intersection number < Tm > for the moduli space of curves 
with one marked point, 

(24^ 

where g is the genus of the curve and m — 3g — 2. We have thus shown 
that the Fourier transform U{s) gives the intersection numbers < > of 
the moduli space of curves with one marked point [5, 18]. 

The replica limit A; — > , where the matrix B is k x k, was studied in 
[6], and it gives the intersection numbers of (12). Note that in the original 
Kontsevich model of (7), the matrix size k was arbitrary, and the universal 
intersection numbers < r„i > are independent of k. 

From other tunings of the external source aj, we may obtain also the 

intersection numbers of the moduli space of p-spin curves [4] with one marked 

point, which exhibit "spin structures". Indeed we may tune the external 

source so that the asymptotic density of states vanishes at an edge as p(A) ~ 
1 

Xp . This will yield the exact values for p-spin curves with genus g and one 
marked point. A spin index j — 0,1, ...,p — 1 is now needed. From this 
tuning of the external source, we obtain the generalized Kontsevich model, 

Z ^ I dBe^^'-'"'^"-'"^^''' (13) 

where B is k x k. The derivation of this partition function from the right 
hand side of the duality formula (5) will be given in the next section. This 
Z has an expansion. 



k. 



m,j,km,j m,j m '^m,j- 
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where 

, I rn— 1 -1 

tmj = i-p)^^^ n (1 + J + IP)^^T^^,- (15) 

1=0 -'^ 

The normahzation constant C in (8) is fixed by (15). The intersection num- 
bers of moduh space of p-spin curves are defined by the integral formula of 
compactified the moduh space Mp_„ [8] 

< TnAUn) ■ --rnMs) >= /- Ct{u) l[ic^iC,)r (16) 

[k + 2)9 JMg,s 

where Uj is an operator for the primary matter field (tachyon), related to 
top Chern class Ct{j^) , and r„ is a gravitational operator, related to the first 
Chern class ci of the line bundle jCi at the ith-marked point. We denote 
Tn{Uj) by r„j, and j represents the spin index (j=0,...,p-l). The indices rij 
and ji are related to genus g and numbers of marked points s through 

(p + 1)(2<7 - 2 + s) = jZimr + 3r + 1). (17) 

i=l 

This intersection theory for spin-curves [4] is known to be related to the 
minimal N = 2 superconformal field theory of Lie algebra Ap^i type, which 
is equivalent to SU{2)f,/U{l) Wess-Zumino-Witten model. A; is a number of 
levels and it is related to phj k = p — 2. This relation is derived from a super 
potential W for the chiral ring structures of primary fields ; we will obtain this 
chiral structure later by the consideration of the n-point correlation functions 
U{si, Sn) (see appendix C). 

As remarked by Witten [4], the limit p ^ —1 {k ^ —3) corresponds to 
the top Chern class without gravity decendants Ci(£j), and this top Chern 
class becomes the orbifold Euler characteristic class [11, 3]. 

For 3 < p, we have to consider the above spin structures for the intersec- 
tion numbers. We find the intersection number with one marked point for 
arbitrary genus, < t„j > [7] for p=3 as 

1 r(2±i) 

^ ^""'^ (12)^f(¥) ^^^^ 

where n = {8g — 5 — j)/3 and j = for g = 3m + 1, j = 1 for g = 3m 
(m-1,2,...). 

In the replica limit. A; — > for the matrix B, a closed expression for 
Uo{si, Sk) is known [6] (the surfix "0" refers to zero external source in the 
Gaussian probability distribution). 
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where a = Yh=i Si, in which Uo{si, s„) is defined by 



Uo{si,...,Sn) = < tre'i-^- • -tre""^ > 

n n k 

= E<ntri?''>n7T- (20) 

li i=l i=l 

Prom (19), we obtain the intersection numbers of p-spin curves with one 
marked point. The rephca hmit A; — > selects only the one marked point 
ribbon graphs on a genus g-Riemann surface. This method gives the in- 
tersection numbers of p-spin curves with one marked point, and the results 
coincide with those obtained from U{s) [6]. 

For two marked points, one deals with the dual quantity U{si, S2), again 
at a critical edge point. The correspondence is the same as for the one- marked 
point. We have 

f/(si,s2) =< tre^i^W^^^ > 
= J dXidX2e''^'+''^^ < tr5(Ai - M)tr6{X2 - M) > 

= lim / d\id\2e'^^^+'^^^——- < [det(Ai - M)]*^i[det(A2 - M)]*^^ > 

= lim / dXidX2e'^^'+'^^^——- < [det(l - iB)f >a (21) 
kiM^oJ 8X18X2 

where A = diag{Xi, ■ ■ ■ , Xi, X, ■ ■ ■ , X2), Xi and A2 are degenerate ki and k2 
times respectively. The matrix B is an {ki + k2) x {ki + k2) Hermitian 
matrix. In the limit of zero replica, one selects the terms of order kik2 in 
the Airy matrix model, for instance a term like -ri-fl) ^-iid we obtain the 
two marked points contribution for the intersection numbers. The Fourier 
transform f/(si, S2), with respect to Ai and A2, gives the intersection numbers 
as coefficients of the Taylor expansion in Si and S2. For the case p=2, this 
was checked for arbitrary genus [5] and it does yield the known values . 

For higher marked points, the argument is similar. For n marked points, 
one considers the terms 

< Tmi • • • 7"m„ > x2mi+l ■v2m2+l . . . \2mn+l ^^^-^ 
^1 ^2 

emerging from the 5-matrix integral. The Fourier transform of these quan- 
tities is given by t/(si, S2, Sn)- 

Thus we have indeed a method for computing the intersection numbers 
of the moduli of curves from random matrix theory, based on the expression 

for [/(Si, ...,Sn)- 



6 



An exact and useful integral representation for U{si, Sn) is known in 
the presence of an external matrix source A [19]. 



(/(si,---,s„) = ;^(tre""...tre-^> 

= jfll^'"- fim-^) de. (23) 

We will use this formula in the following sections to obtain the intersection 
numbers of p-spin curves of arbitrary genus g for n-marked points, through 
an appropriate tuning of the external source aj, in a scaling large N limit. 

3 The p-dependence of the intersection num- 
bers with one marked point 

The partition function Zp for the generalized Kontsevich model is given by 
the k X k Hermitian matrix B, 

Zp^l dSe^*"-^'^'-*--^^ (24) 

This model is obtained, after use of the duality, from the expectation val- 
ues of characteristic polynomials (5). We take here an external source A 
with (p — 1) distinct eigenvalues, each of them being times degenerate : 
A = diag(ai, ...,ai, ....,ap=i, ....,ap_i). After duality, the expectation values 
of characteristic polynomials become 

P—I P— 1 ■ p— 1 

< n det(a,-^5)5^ >=< expl^^ trlog(l-— )+7VX^log(n > (25) 

a=l a=l a=l 

We now specify the (p — 1) distinct eigenvalues of the external source by the 
(p — 1) conditions : 

E;^ = 0, (m = 3,4,...,p) 

(26) 

Then, the expectation values of the characteristic polynomials lead to (24) 
in the double scaling limit. 

We first consider the intersection numbers with one marked point. They 
are related to the coefficients of tr-^, in the zero- replica limit A; ^ as we 



P-1 I 

a=l 
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have seen in the previous section. In this hmit the matrix A can be taken as 
multiple of the identity A = A ■ 1. We introduce a coupling constant g as 

B^- (27) 

9 

Z,^j dBeT^'-'"'"'--^'''' (28) 

We set Qs = g^^^ and tune A so that S's = f • Then we obtain, after the shift 
B^l + B, 

Z,^J dBeT^''^'^-^^'^'-^''^'^''^ (29) 
Expanding for small B, we have 

/P 4-p2 I P(P-l) <-,p3 I P(p-l)(p-2) <.,p4 . . 

dBe^^"^ +^Ti^*'^^ + — 5ii — ••• (30) 

We now expand in powers of Qs after the replacement B i^J^B. 

Using the replica formula (19) for one marked point, we obtain limfe_+o j: < 
trB^ >= l,limk^o ^ < (trB^)^ >= 3, .... Using these values for the products 
of vertices, we have 

k^o ^ P 24 >^ 1920 p3 ^ y^sJ 

(31) 

The coefficients of the above expansion are intersection numbers multi- 
plied by tmj in (15), the first term for genus one, the second for genus two, 
etc.. From (15), we have ti o = — ^tr^^^. Therefore, the intersection number 
of one marked point for genus one becomes 

< Ti,o (32) 

For genus two, we have ^3^2 = (— p)~^3(3 -l-p)(3 + 2p)tYj^^, and 

^ (p-l)(p-3)(l + 2p) 
< >^=^^ p. 5!. 4^. 3 ■ ^^^^ 

For g — 2 and p — 2, we have ^4,0 — ~^^tr^ and this gives < t^^ >g=2— 
■jj^. The expansion of (31) can be obtained for any higher order of genus g 
by the use of replica formula of (19) although the evaluation becomes tedious. 

We now turn to the dual model, formulated with N x N random matrices 
M ; the Fourier transform U (s) of the one point correlation function is given 
in (23). 

We are still in the case in which the external source Ua takes values p — 1 
different values of ai, 02, Op-i with degeneracy. 
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We have from (23), 



„2 



Uis) = — / — e-eSr=i ^ iog(a.-«-.)-Eri ^ losK-) (34) 
^ ' NsJ 2m ^ ' 

Expanding u{s) for small u and s, we have 
^^^^ " TVs)' 2i7r 



With the conditions (26), we obtain, after the shift u ^ u — ^, 



(35) 



U(s) = —[ i!ie-^[("+^^)^^'-("-^^)''^'] (36) 
^ ' NsJ 2m ^ ' 

where c=^E^- 

Expanding the exponent, we obtain 

This integrals yield Gamma functions after the replacement u— (— )^''^, 



1 1 /""^ 1 1 

U{s) = —-^ / dtt-^-\-' 

^ ' Nsp-K (csy/P Jo 



Nspn (csy 

p(p-l) 2+§ p 1-| 2)(p-3) 4+1 i 1-^ 

xe 314 * 5532 * + 

(p-l)(p-3)(l + 2p) _ 3 
5! -42 -3 ^ ^ 



(P - 5)(p - 1)(1 + 2p)(8p2 - I3p - 13) 3 



rr(i - 3) 



+ {P - 7)(p - 1)(1 + 2p)(72p^ - 298p2 - 17^^ + 562p + 281) 
with y = cps^"*"?. 

Comparing this expansion with (31), we find a common p-dendence, but 
(31) has additional factors. These additional factors should be included in 
the normalization. The intersection numbers < t„j >g are the coefficients of 
(38), and then the two results of (31) and (38) coincide. The expansion of 
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U{s) is easily obtained up to arbitrary order of genus g since U{s) is simply 
given by (36). Thus the dual model is simpler than the the partition function 
of the matrix B. 

The expansion of U (s) in (38) is genus expansion. We write this expansion 

as 

U{s) = E < rn, >. ^r(l - 1 - ^)c^p-^.(^-^)(^+|), (39) 

where n and j are given by 

{p+l){2g-l)^pn + j + l. (40) 

Prom (38), the intersection numbers < Tnj > are determined explicitly for 
arbitrary p. 

For g — I case, by the condition of (40), we obtain n — 1 and j — 0, and 

p-1 

< Tl,0 >g=l^ (41) 

For g = 2, we have n — 3 + The intersection numbers for arbitrary 
p become 

(p- i)(p-3)(i + 2p) r(i-^) 

p-51-42.3 r(l-i±i)' 



< TnJ >g=2- „ ,2 O 7^71 T+j^' l^^J 



For instance, we obtain < r3,2 >g=2— ^^'"^ p5U2.3 "^^^^ ; which agrees with the 

result of (33) evaluated from (31). For p — 2, we have from this formula, 

< ^4,0 >g=2= (24p2!) siuce the ratio of the gamma functions becomes —2. 

For 5f = 3, we have the intersection numbers for arbitrary p with n — 

5 + ^, 
p ' 

_ {p-5){p- 1)(1 + 2p){8p' - 13p - 13) r(l - I) 

^ ^9=3- p2 . 7! . 43 . 32 _ l±i) ■ V^"^; 



For 51 = 4, we have the intersection numbers with thecondition n = 
{p -7){p- 1)(1 + 2p)(72p^ - 298p3 - ITp^ + 562p + 281) 



7+^, 
p ' 



p3 . 9! . 44 . 15 



Ffl 



X 7^7^^- (44) 



r(i-^) 



Thus the expansion oiU{s) gives the intersection numbers of one marked 
point for arbitrary p in the case of fixed g. We have given the explicit 
expressions up to order g — A only. Also the integral representation of U (s) 
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in (36) gives the intersection numbers for arbitrary g for fixed p. As sliown 
before in [7], the intersection numbers for p = 2 becomes 

< ^"'0 (^^ ' = % - 2) ■ (45) 

For p = 3, 

1 r(i±^) 

^ (12)9^! r(4^)' ^^^^ 

with 3n = 8g — 5 — j. 

We now consider the interesting hmit p — > — 1 . We consider the spin 
index j as j — in this hmit. Prom the condition of (40), we obtain n — 1. 
The intersection number of p = —1 case is then < ti,o >g, which we write 
simply as < T >g in the following. Prom the previous evaluations in (41)- (44), 
we obtain in the limit p ^ — 1, 



< r >a=2= 



p- 1 _ 1 

^4 ' "12 

(p-i)(p-3)(i + 2p)r(i-|; 



5!42 r(l - i) 120 

<"^-^^-2i0- 
These numbers are the Euler characteristics x{Mg,i) [3, H]- 

x{Mg,,) = ai-2g) = -^ (48) 

where C is the Riemann zeta-function and B2g is the Bernoulli number; B2 — 

1 R — 1 R _ 1 R _ 1 

6' ^4 - -30, ^6 - 42, -«8 - 30 ■ ■ •• 

The logarithmic term of the Penner model follows indeed in the limit 
p ^ —1 of (36). Then computing U{s) for p — > —1, one obtains the Euler 
characteristics. 

In the limit p — > —1, c = ^ J2 — fft is N, and from (36), C/(s) is given 

by 

U{s) = — / — e "-5" 



Setting 



TVs 7 2i7r 



-y ^ l + e y 

= e ^, (ti = ) (50) 

u+1 ^ 1-6-2^^ ^ ^ 
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one has 



N J 2n {1 - e-yy"" ~ Jo 2nl-e-y ^^^^ 

Noting that 

1 oo j.n—1 

E V (52) 



we obtain U (s) as 

- ^ ^ 1111 

~ ^ ~ 27V ^ 12A^2 ~ 1207V4 ^ 252 7V6 ^ " ^ 

Then we obtain the genus g, orbifold Euler characteristics x{Mg^i) — C{1 — 
2g) — —^B2g from the term of order Thus the analytic continuation 

for negative p holds for the dual model. 



4 The n-point correlation functions 

We consider the two-point correlation function U{si,S2) defined in (1). Not- 
ing that the two terms of the determinant in (23) become, after the shift 



1 1 



U1-U2 + ^{Si + S2) U1-U2- ^{si + S2) 

= ( V \-r n)^^ (54) 

Ul - U2 - i^[Si + S2) U1-U2 + ^{Si + S2) S1 + S2 

we write it as 

1 1 



U1-U2 + 2^(S1 + S2) U1-U2 - 2^(Si S2) 

= — / dxe--("i-"^)sh( — (si + S2)) (55) 
Si + S2 JO iJS 

We have at the same p-th critical point defined for the one point function, 



2N 1 f°° f 1 

U{si,S2) = ■ ,^ / dx / duidu2sh{—x{si + S2)) 

Si -|- S2 {Imy Jo J Zis 

xexph^ E ^(E(«^ + ^s^r^ - Y.{u. - ^s.Y^-^) (56) 
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Wc use the notation, c = After the change of variables Ui ivi 

(i=l,2), the rescalings Vi and x we obtain 

2N' ,1,1/"°°, dvidv2 , , X i , 
si + so So Jo J-oo (27rr 2iV' 



—ixvi+ixv2( — )P 

e "2' 



1 2 



n G{v,) (57) 

i=l 

where we used A^' = N{^^)p, and [|] = | for even p and [|] = for odd 
p. The factor G{vi) is given by 

G(^;j-exp[ ' (2m + l)!22-(p - 2m)!7V'2m^^ J-^^^^ 

The genus g of the terms in the expansion (58) is given by the exponent 

of ^lyja?- We are interested in the terms of type ^ S2 ^ in (57). The 
correspondence with the variable t 

n,m ~ \pn+m+i is 

~ t„,m (59) 
Thus we obtain the intersection numbers < Tnj^,miTn2,m2 >g from the coeffi- 



m- 



-inill n2+^^i2±i _ 1 2+i 



cients of the terms ^2 ^ ■ instance, we have from s[s2 

(57), 

p-1 

< To,oT2,0 >g=l= (60) 

This value coincides with < ti^ >g=i= and we have 

< To,oT2,0 >g=l = < Ti,o >g=l (61) 

which is consistent with the string equation. Indeed the generating function 
F for the intersection numbers satisfies the string equation [4], 

dF 1 , °° dF 

W. = ^ X! ^'""'^O.m^O.m' + X! X! Wl,m^T (62) 

t/ro,0 ^ m,m'=0 n=l m=0 ^'^n,m 

where the metric rf^'^' = 5m+m',p-2- 

If we substitute F = ati^Q (a is some constant) in the second term, then 
we have = ^2,0^^^ = at2,0) which yields ato,o^2,o after integration. Thus 
we have < to,oT2,o >=< ti,o > from the string equation (62). 

Note that in the two point correlation (two marked points), there is no 
genus zero contribution to the intersection numbers, an easy consequence of 
(57). 



13 



It is useful to define the higher Airy functions 4>p{x) by 

M-) = J |e-f---^ (63) 
When p = 3, it reduces to the usual Airy function Ai{x), 

(j)3{x) = Ai{x) = — / cive3^ +*^'' (64) 

ZTT J-oo 

which satisfies the differential equation, 

U^y = xu^) (65) 

For general p, we have 

j,^p-i = ^M^) (66) 

Expanding in powers of after integration by parts for the function (f)p{x), 
we obtain all intersection numbers from (57) for two marked points. In the 
case p = 2, the function (f)p{x) is Gaussian , and we are led to 



fi + «2 ^„m!(2m+l)^ STV^ ' 



which has been obtained in [5, 18]. The intersection numbers are given by 

N2g 



U{Si,S2)^ J2 < Tni,0rn2,0 >5 (68) 

ni,n2 



where the genus g is specified by 

m + 712 = 3^ - 2. (69) 

For the p = 3 case, we have 

C/(si,S2) = -— -(^)^ / dxsH--s!{si + S2))A{x)A{-x{^)-3){70) 

Si + S2 yS2 JO z.l\ S2 

which gives the intersection numbers of two-marked points. For instance, we 
obtain ^ ^ 

< To,oT2,0 >3=1= "^I'O >fl=l= ^^^^ 

In general, the integral formula of (23) gives the n-point correlation func- 
tion U{si, Sn)- In Appendix A, we compute the intersection numbers for 
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three marked points from the three point functionf/(si, S2, Sa)- In Appendix 
B, the intersection numbers with four marked points are computed from the 
four point correlation function. 

The intersection numbers for the primary field Uj of (16), in the genus 
zero case, are particularly important, since they have algebraic structures 
related to a super-conformal field theory [4, 13]. They are expressed as 
< U.m=l ^o,qm >• 

In Appendix A, we compute the intersection numbers (133) as 

< 70,91 7"0,q2 "^0,93 >5=0= <^9i +92+93,^-2 (72) 

Prom this result, the free energy F follows 

P = ^< 'n),qiT0,g2 70,93 >g=Q ^0,91^0,92 ^0,93 + O(t^) (73) 

To make the algebraic structure more explicit, we define the structure con- 
stants Cijk by 

" Om^dtk ^^^^ 

where 

= to,i_i, (i = 1, - 1) (75) 
For instance for p=5, we have 

F = ^^0,0^0,3 + to,oto,ito,2 + + 0{t^) (76) 

and the striicture constants are C114 = C123 = C222 ~ 1- Prom these struc- 
ture constants, following [4] one can construct the super potential W (see 
appendix C). The free energy F, which generates the intersection numbers 
of the moduh space of p-spin curves, has been conjectured by Witten to be 
solution of the Gelfand-Dikii hierarchy [4] . We present in Appendix C, this 
Gclfand-Dikii hierarchic equations as well as the construction of the super 
potential for the primary fields. 

Thus, we find that the intersection numbers, derived from the integral 
representation of [/(si, s„), satisfy indeed the Gelfand-Dikii equations. 

Returning to Witten's conjecture, we note that the definition of the in- 
tersection numbers by the vector bundle integration over the compactified 
moduli space Mg,n in (16), is similar in structure to the integral representa- 
tion for t/(si, Sn) of (23), although t/(si, s„) involves a summation over 
all genuses for fixed n-marked points. 
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5 Time-dependent Gaussian matrix model 



Let us briefly recaU how one shows that the time-dependent one matrix model 

is equivalent to a time-independent two-matrix model when the distributions 
are Gaussian. ( In an older work we had considered the time-dependent 
Gaussian matrix problem, and computed time-dependent correlation func- 
tions [20]). 

The time-dependent Gaussian matrix model is the partition function (i.e. 
the Euclidean path integral) for the matrix quantum mechanics, with action 



+ M^) (77) 



where M isa,n NxN Hermitian matrix ; (the dot stands for time derivative) . 
This model, at criticality, is known to describe gravity coupled to matter of 
central charge c = 1 [21]. 

The time-dependent correlation function is defined by 

p(A, t) =< ^tr5(A - M(ti))^tr5(/. - M{t2)) > (78) 

and, from time-translation invariance, is function oit — 1^2 — | ■ The Fourier 
transform of this quantity is 

U(a, /?) = 4^ < tre^^A^^^^^tre*'^^^*^) > (79) 

This correlation function may easily be reduced to the correlation function 
of the time-independent two matrix model in the Gaussian ensemble [20]. 
Indeed this harmonic oscillator quantum mechanics leads to 



Rescaling of A,B,a and /3 by a factor Ve *sinht, we obtain a time-independent 
two-matrix model 

U{a,(3) = ^ / (i>ldStre^"^tre^^^e-^*'-(^'+^'-2^^^), (81) 

with a coupling constant 

c = e-K 

For convenience, we denote A and B by Mi and M2 in the following . 
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6 Duality formula for the two-matrix model 



We consider the correlation function of the characteristic polynomials in the 
two-matrix model, 

J =< n det(A<, - Ml) n detiiip - M^) > (82) 

a=l /3=1 

where the average is perfomed over a two-matrix Gaussian distribution with 
an external source A acting on one of the two N x N matrices, 

P(Mi,M2) = le-^t'-^l'-5t'-^2-ctrMiM2-trMiA ^gg^ 

The external source A will be used here again mean to tune a p-spin structure 
in the moduli space. Note that when i — > oo the parameter c = e"* vanishes 
and the two matrices Mi and M2 decouple. 

From (82) we shall determine a new dual model of Kontsevich type in the 
large N limit. The duality formula for J is obtained by the use of Grassmann 
variables as in the one matrix case [6, 16] , but for the two matrix case a new 
structure appears. 

Let us introduce the Grassmann variables ipf and xf , where a — 1, ki, 
and (3 — 1, k2. Then 

J =< / <i^,V'**e"l*"<^->*-«»<«-"='«> > (84) 

Since the probability P is Gaussian, one can integrate out the matrices Mi 
and M2. This generates four-fermion terms that may be disentangled with 
the help of three auxiliary matrices : Bi a ki x ki Hermitian matrix, B2 a 
k2 X k2 Hermitian matrix and D a complex ki x k2 rectangular matrix. The 
identities 

= I d5,e--"^^^7n?"^^^^ (85) 

allow to represent J as 

J = j rf5irf52ci/^^rf/^e-Tt^(^!+^i+2Dti?) 



N 

X JJdet 

i=l 



{^a - iZ^)(^a,a' + 7T=f "^1 \/ 



(88) 
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After the shift Bi -Bi+^Vl — c'^^a,a'Sa,a' and B2 — >■ -B2+^\/l — c^lJ'i3,i3'Si3,i3', 
we obtain the dual expression for J : 

X e-Stit'-i°g(i-^*), (89) 
where the matrices Xi are defined by 

/ iv/r^^fi /c(i-c2)^ \ 

^ c c ai ' 

Wc now expand log(l — Xj) in powers of Xi. Let us first consider the case 

of a source matrix A mutiple of the identity : = a, i = 1 ■ ■ ■ A^. Imposing 
the constraint 

a = (1 - c') (91) 

the quadratic term trSf cancels with the one coming from the expansion 
of log(l — X). This critical constraint corresponds to the edge of the the 
spectrum for the matrix Mi. Note that the term is not cancelled at this 
critical point because of the coupling c. 

Given the factor N in the exponent, the edge scaling limit under consid- 
eration corresponds to 

El - 0(X^^), B2 - 0{N-^), D - 0(iV-5) (92) 

in the large N limit. In this limit most terms disappear ; for instance 

Ntr{D^DB2) ~ A^"^ (93) 

is negligible. Then, in the large N limit (92), we obtain the partition function 
Z, i.e. J after dropping the negligible terms, 

(94) 

Since the matrix B2 matrix is decoupled we can integrate it out. Then, 
dropping the decoupled part, we find the partition function 

J dBidDUDe-'''^'^'+-3''^'+'*'^^^^' (95) 

where we have absorbed the powers of N given by the scaling (92). We may 
now integrate out the matrices D and ; this yields a one matrix integral 
with a logarithmic potential, 

Z^ J dBie-^'"^'-'''''^''^^'-'''^'^'. (96) 
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The appearance of a logarithmic term is a characteristic of models with 
central charge equal to one. 

We now consider the free energy of this logarithmic Kontsevich model 
(p=2) (96). Three different, but consistent, methods will be used. (For 
convenience k2 is denoted as q in what follows.) 

i) HarishChandra-Itzykson-Zuber method 

After use of the HarishChandra-Itzykson-Zuber formula, the partition 
function Z is given by 

Z ^ j (iBe^^''^^~i^'^^°^^~'^'^''^^^ 

^ W)!% '^'"'^^''^ ^ ^i^ii^H^-<^^i (97) 

where the XiS are the eigenvalues of B, the U the eigenvalues of A, and A (a;) 
the Vandermonde determinant A(a;) = Y{i^j{xi — Xj). It may be replaced in 

d 



the integrand by the Vandermonde of differential operators ^ 



n(g^(|-^)nc(y (98) 



where 



Cik) = J dxet^'-^^'?-«i°s^ (99) 
Rescaling Xi — > Xi/2^/^ and k and with the change I — > il, we have 

0) [ rfa;e-^"'+iI^"'-^^°^^'+^) (100) 



11+ 

Expanding for large I, we obtain 

log^ = -[ltl + ^t3 + qt2ti + ^qH3] 

+[^*?^3 + + ^^1^5 + ^qh + qt\U 
13 1 

+qtit2h + -^qtl + -qHih + -^q^tl + qH2h 

+\q%\ + o{^) (101) 

where we have used the moduli parameters 

fei I 



1=1 
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When q ^ 0, we recover the resuh of the one-matrix Kontsevich model. 

For the relation to the genus g, we have to identify the powers of In 
the limit in which 

N, k ~ A^3 (103) 

we find 

~ 0{Nl),tn ~ 0(A^^~5") (104) 
The genus expansion of the free energy 

oo 

\ogZ = Y.agN^-^^ (105) 

3=0 

follows from this limit. For instance, tl,qtit2, and qH^ are contributions to 
genus zero, and to genus one. 

ii) replica method 

We return to the integral (97) . After the shift B B + A, which 
eliminates the terms linear in B, one can expand for large A. Then the 
logarithmic potential, trlog(B + A) expanded in powers of A^^ yields trB'^ 
vertices. The situation is similar to that of the generalized Kontsevich model 
where the trBP"*"^ terms led to the moduli space of p-th curves, and spin 
structures appeared. The occurence of t2i^4 and indicates this fact. 

We consider the moduli space for Riemann surfaces with marked points. 
Although there is a logarithmic potential, the model allows one to consider 
marked points, whose number is equal to the number of tn- 

We have developed the replica method /ci — >^ in a previous article [6]. 
Any average of the products of vertices trB"^ are obtained in the replica limit 
ki 0, where B is a ki x ki Hermitian matrix ; for a Gaussian ensemble, 
the average is given by the replica limit formula (19). 

Expanding the logarithmic term, after the shift B ^ B + A, and the 
use of the formula (19), we obtain the replica limit, which gives the required 
intersection numbers with one marked point. Indeed it leads easily to 

log^ = + ^^')^3 + + lq')te + O(^) (106) 

This result agrees completely with the expression (101) for one marked point. 

iii) differential equation of Virasoro type 

Since the free energy logZ is expressed in terms of the moduli parameters 
tn, as was the case in the original Kontsevich model (q=0 case), it is natural to 
investigate here again the KdV-like differential equations or string equations. 
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Although there is a logarithmic potential one may use a Schwinger-Dyson 
equation [22, 23]. 

We first consider the simple case, ki — 1, a one by one matrix, i.e. a 
c-number. Denoting li — x, one finds 

Z^e-'^^'^g{x). (107) 

The Schwinger-Dyson (Virasoro) equation follows from the identity 

I ^5^gttrB3-itrBA2-gtrlogB _ q ^^gg) 

The matrix B is the replaced by For the logarithmic potential, it 
means (j^)^^- Therefore we need to apply a differential operator in order 
to get rid of this integral. We find easily that when B is just a real number 
{ki = 1), the function g{x) satisfies a third order differential equation, 

[(J_,3_2 + 2,-.Aj£Z^,M.o (109) 

i.e. 

(-(1 + 2g)(5 + 2g)(9 + 2q) + (-10 - 48g - 2A(f )x^)g 

+ ((66 + 96q + 2V),T + (24 + A8q)x^ + 16x^)^'(a;) 

+ (-36a;^ - 24:qx^ - 24a;^)/ + 8x^g'"{x) = 0. (110) 

This provides the large x expansion, 

/ N 1/5 q^. 

1 , 385 73 161 2 7 o 1 

- ^^^1152 + 24^+ 4^^ ^e'^'+S^^ 

1 ,85085 6259 58057 ^ 2075 o 

+ — f \ a -\ H 

82944 384 ^ 2304 ^ 144^ 

For /ci = 2, a two by two matrix, we denote x — li and y — I2 ; then 

xax xox x^ — xox yoy ox (xy^^^ix + y) 

(112) 

The solution, after symmetrization over x and y, agrees with the expression 
(101). 
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When X is small, the equation (110) leads to a different series expansion. 
The differential equation for g in (110) has three different solutions, 

g{x) ~ x'^+^,g{x) ~ x''+^,g{x) ~ (113) 

For small x, we have from the first solution, noting that I — x, 

Z = e-¥\l + 1/3 + (__L^ + J_ye ^ 0(^9) (114) 

This solution for small Z is in a different phase from Kontsevich's phase ; 
it may be related to one of the two phases of the unitary matrix model 
[23, 24, 25]. 

We now consider the case p > 2. After the integration over the In- 
fields within the matrix X, we also obtain the logarithmic term trlog-B, but 
corrections appear. By tuning the external source with the conditions (26), 
we obtain 

Z^l dXe-^*'-^^^'+*'-^^^ (115) 

where X is given by 

x=(^, ;^). (116) 

where we have scaled out the factors \/l — in X^ of (89), and put B — Bi 
and ^2 = 0. We expand the potential, 

tTXP+^ = tTBP+^ + {p+ l)tTDD''BP-^ + l(p + l)(p - 2)tT{DD'')'^BP-^ 

+ ^(P + 1)(P - 3)(p - 4:)tr{DD^fBP-^ + ■■■. (117) 

For p = 3, we obtain, 

Z^ J (/5^i^t^^e-t^^^'+*'-^^'^'+^*'-(^^')']+*'-^^'. (118) 

The integration of the D-field can not be done explicitly for the general 
values of ki and A;2 {B is a ki x ki Hermitian matrix and D is a fci x A;2 
complex matrix). We make here a perturbation for the large B in lower 
orders. Expanding the term exp{—^tY{DD^y), we find 

For p — 4 case, the partition function Z becomes similarly 

= y ^^gtrBA4-itrB5-3fe2trlogiJ-fc2tr;J^-fc2tr^tr^+0(;g^)_ (^20) 
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For general p, after integration over the D-fie\d in a perturbation, we 
obtain 



In order to identify the intersection numbers, we expand in powers of j^. 
For this purpose, we shift B ^ B + A. New terms in the exponent, which 
are corrections to the logarithmic term, have the form of the product of two 
traces. In the large A case, the fifth term is of the form 

tryi T^^^T-. 7^r^ 2(tr— S) • tr- — - H (122) 

The term tr^^^ is tp_i ; such terms appear in the c = 1 string theory [14, 15]. 

We now evaluate the intersection numbers for a small number of marked 
points, and for lower orders in In this case, we can neglect the above 
correction terms, and we approximate the partition function by 

Z = J (^5e-^trBP+i-(p-l)fe2trlogiJ+trBAf ^^23) 

We find, at order j^+i, up to three marked points, 

ICS. - 



+ {higher order) (124) 

where the overall factor ^ is a normalization constant absorbed in A. Thus 

p 

wc find the two-matrix model for c = 1, obtained from the characteristic 
polynomials, reduces to the one matrix model, and the topological invariants 
becomes similar to the intersection numbers of p-spin curves. 

When we put p — > —1, we find that the last term in (121) (tr-^)(tr-g^) 
behaves like 

= - (125) 

where we make a shift S — > A + S. In the case p — > — 1, when B is order of 
A, all the terms of the potential should be order of one, and indeed tm(tri?"^) 
is order of one for the A-dependence. Therefore, the potential has a series of 
^m(tr-B™). Such terms appear in the c=l string theory from the calculation 
of the tachyon correlators [14, 15]. 
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7 Discussion 



In this article, we have considered the exphcit p-dependence of the inter- 
section numbers of moduh spaces of p-th spin curves based on one and two 
Gaussian matrix models. 

In the one-matrix case, the limit p — > — 1 gives a generating function 
of the Euler characteristics. In the two-matrix case, we have obtained a 
logarithmic matrix model with polynomial corrections, which is related to 
the generating function for the tachyon correlators [14, 15]. 

The duality, on which the present analysis relies, is the relation between 
the characteristic polynomials of two different Gaussian matrices. The char- 
acteristic polynomials are computed as determinants, expressed in terms of 
Grassmann variables = l,...,N,a = l,...k). In the large N limit the 

p-th singularity is tuned through an appropriate choice of the eigenvalues 
of an external source matrix A. One parameter remains, namely the num- 
ber of different A^, a — 1- ■ - k. The Fourier transform with respect to the 
Xa yields the correlation function C/(si, s„). The symmetry between A'^ 
and k becomes then implicit. Although this duality might be related to the 
open/closed string duality [26, 27, 28], we have not been able yet to reach a 
clear picuture in this direction. 
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Appendix A: Three point correlation function U{si, 82,83) 

For the three-point correlation function U (si, 82, 83), we address ourselves 
to the determinant terms in (23) similar to the two-point case. The longest 
cycle in the determinant of a 3 x 3 matrix is 



det{aij)\iongest = 012023031 + 013021032 (126) 



where a,,- — 



1 



^-^ Ui-Uj + ^{Si+Sj) ' 

We consider the first cycle of (126), ( the second cycle is almost the same). 



U1-U2 + |(S1 + 82) U2-U3 + \{S2 + S3) U3-U1 + ^(S3 + Si) 
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POO poo rj^ 

/ dx dysh{-{si + S2 + S3)) 
Jo Jo z 



Si + S2 + S3 

We express the two terms as 

U(s,,S2,S3)^U' + U''. (128) 
After the shift Si ^ j^, using the notation A^' = N{^0)p, we have 

11 1 I 



U'' = ^ ^ (-)^ / / C^?/sh(-^«r(.l + .2 + «3)) 

Si + S2 + S3 S3 h Jo 2N' 

11 11 

where G'(i'i) is defined by (58). Expanding G{vi) in powers of U{si, S2, S3) 
is expressed in terms of the function (j)p{x). 

The intersection numbers < Tn^,miTn2,m2^n3,m3 > is obtained from the 

771-3 + 1 

coefficients of sT ^ 3^^""'+^^ s^^ ^ . 

In this three point correlation function, non-trivial genus zero terms ap- 

pear. From U^^ in (130), we obtain the term sf sf S3 ^ in the large A^' limit. 
This leads to 

< T-0,0T0,0T"0,p-2 ><y=0= 1 (131) 

Since there is terms of {^)py and {^)px in (130), these terms contribute 

l+gl 1+92 1+91+92 

in the large N limit as ^ S2 ^ S3 ^ , and we obtain the intersection 
numbers, 

< TO,giTo,q2ro,p-2-gi-g2 >g=0= 1 (132) 

This is related to the property of ring correlators found in [4] 

< To,qiTo,q2To,g3 >g=0= <^qi+g2+93,p-2 (133) 

which is important for the chiral ring theory and superconformal theory for 
the primary fields. Prom this result, the generating function F is obtained 
as 

= X] < 'n),qiT0,g2'^0,g3 >9=0 ^0,91^0,92^0,93 + <^(^^) (134) 
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and the superpotential W can be constructed from the structure constants 
Cijk defined by 

= dUdtjdt, ^^^^^ 

where we put 

ti^to^i-i, (i = 1, - 1) (136) 

If we consider only primary field, neglecting gravitational descendants, 
we only need the terms Ilm^o.m- When we consider this primary field, in the 
genus zero case, we obtain therefore ,for instance for p=5, 

F = ^^0,0^0,3 + ^0,0^0,1^0,2 + ^^^,1 + 0{t^) (137) 
and the structure constants become Ciu — C123 = C222 — 1 for p=5 case. 

Appendix B: The n-point correlation function for n> A 

The calculation of the n-point correlation function C/(si, s„) at edge 
singularities follows the same steps as for the n—2 and 3 cases. For the 
discussion of the higher chiral ring structure , we need more than three 
points, and thus we consider n>4. 

One of the longest cycle terms in the determinant for the four point 
correlation function U (si, S2, S3, S4) is 

012023034041 = ■ 7 ■ (012 + 023)(034 + O41) 

•Sl + -52 + S3 + Sa 

1 1 ^ . ^ 

X w 7, ^ T7 7, ^) 138 

U1-U3 + + 2S2 + S3) U1-U3 - ^(Si + 2S4 + S3) 

with ^ 

Ui-Uj + \{Si + Sj) ^^^^^ 



This term can be expressed by the integrals 

2 

Sl + S2 + S3 + S4 



012023034041 = ^ /°°dxdyd^e-^("i-"3)-^(^2-^4)x 

Sl + S2 + S3 + S4 ^0 

xsinh(-a;(si + S2 + S3 + S4)) 

x[exp(-^y(si + S2) - ^2;(s3 + S4) - uxy + yui - zu^ + zua) 
+exp(-iy(si + S2) - ^^(si + S4) - uiy + yu2 - zu^ + zui) 
+exp(-^y(s2 + S3) - ^2;(s3 + S4) - U2y + yuz - zus + zu^) 
+exp(-^y(s2 + S3) - ^z{si + S4) - U2y + yus - zu^ + ^-ui)] (140) 
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Using the same change of variables and scahngs as before, we obtain 

U{si, S2, S3, S4) ^U' + U'' + U''' + U'^ (141) 
These four terms are given by cr = Si + S2 + S3 + S4, 

1 i 1 i 1 i 

xexp[-^(s2 - S4)sf X - j^{si + S2)s^y - ^(sa + Si)slz 
—ixvi — i{—)pyvi + iyv2 + i(— )pa;t'3 — izv^ + i(— )p2;t'4) (142) 

Sl S3 S4 



1 i 1 i 1 i 

xexp[-^(s2 - S4)s[a; - ^(si + S2)s^y - ^(^i + S4)s|z 
— ixvi — i(— )pyi'i + i{—)pzvi + iyv2 + i{—)pxv3 — i{—)pzv4) 

Si Si S3 S4 



(143) 



^ = -^^^>y (2^""^^w^^"\n^^(-o 

1 i 1 i 1 i 

xexp[-^(s2 - S4)s[a; - ^(^2 + S3)s|2/ - ^(^3 + S4)s|z 

—ixvi — iyv2 + i{—)pvsx — izv^ + i{—)pyvs + i{—)pzv4) (144) 
S3 S3 S4 



jy 2N'\1 1 f dv, X 1 ' 

1 i 1 i 1 1 

xexp[-^(s2 - S4)sra; - ^(^2 + S3)s|y - ^(^i + Siis^z 

, v-^sa . , ..si 1 S2.1 ..S3 1 

Sl S3 S3 S4 

(145) 

2 2 X— i 

Prom U^^^ , we obtain the term sfs2S3 ^S4 which gives the intersection 
number < To,iTo,iTo,p_2To,p_2 >s=o- In this large N' limit, we have 

jjiii ^ /■ / 7T^^ " (tt^i ^) • (^^ssly) • i{—)pyv3 

S4 Jo J (Znr 2 2 S3 
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IP 53 1 

exp[ '^v^ — ixvi-\- i{ — )pzvi — iyv^^ x 

V i si 

r..Sl ..S3 1 

exp[z( — jxvz + z( — jpyv^ — i\ — )pzv/^\ 
S3 S3 S4 



(146) 



Expanding the factors exp[i(^)a;f3 + i{^)'^yv2, — 'i'{^)'^zvi\ we obtain the 
series of the intersection numbers for the primary fields in the genus zero 
case, 

< 7"0,qi7"0,q2'^0,p-qi-g2+«3> '^0,p-2-q3 >g=0 (147) 

where gi, ^2 = 1, 2, and ^3 = 0, 1, 2, .... The other three terms [/^, U^^ and 

gl+l 92+1 91+92^93 93 + 1 

C/'^^ do not yeild terms of the type Si S2 S3 ^ S4 . 

For p=2, the term (147) does not exist, since ro,i is not aUowed. For 
higher n-point correlations (n > 5), there is no correction for the same reason. 
Therefore, for p=2, the function F for the primary field is 

F=\tl, (p = 2). (148) 
For p=3, we obtain from (147) and (133), 

F = ^4o^o,i + (149) 



For p=4, we obtain (147) and (133), 

F = 2^0'0^o,2 + 2^0.0^0,1 + Y^'^0,1'^0,2 + ^—^^0,2 (150) 

The last term is evaluated from the five point correlation function, which has 
the form, for general p, 

^ 1 + 91 -^+92 I 92+1 a+91 1^1 

^1 ^ =^2 ^ =^3 -^4 ^5 ^ ^ ^0,p-2-(jii^0,2+g2^0,p-2-g2^0,2+gi^0,p-2 (1^1) 

which leads to the last term in (150) for p=4. We have investigated the 
intersection numbers of primary fields, but other gravity descendants can be 
obtained in the same ways, which would be in factor of tn^m (^ 7^ 0). 

Appendix C: Ginzburg-Landau potential for primciry fields and 
Gelfand-Dikii equation 

The structure constant Cijk defined by (135) are obtained from the n- 
point correlation function through the intersection numbers with n marked 
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points. In this appendix, we discuss the relation to the superpotential [4]. 
Using the notation 

ti^to,i-i (i = l,2,...) (152) 
and the metric 77""^ = (^n+TO,p; we define 

C^, = E QjmV"''. (153) 
m=l 

In this notation, F becomes, in the p=4 case for instance, 

F = ^tlt, + U,tl + ^tltl + ^tl (154) 
We find that the Witten, Dijkgraaf, VerIinde,@VerIinde relation® [4, 13] 

ij ^mkl — ^ik^mjl \^^^) 

holds for the structure constants that we have computed. 
Then the C^j have a ring structure, 

Mj = E {mod[W'{x)]) (156) 

k 

where 0j is defined by the derivative of the Landau-Ginzburg potential W{x) 

dW 

We have obtained the function F by the evaluation of the intersection 
numbers of primary fields up to the 6-point correlation function. For the 
p=5 case, 

^ = 2^0'°^°'^ ^0,0^0,1^0,2 + g^o,i + 4^0,1^0,3 

+ 0^0,1^0 2^0,3 + 7:^0,1^0 2 + 0^0 1^0,2^0,3 
/ O Z 



1 ^ 1 2 3 -'-6 
+ ^2^0.2 + g ^0,2^0,3 + Y20^°'^ i^^^) 



This leads to C213 = C411 = C222 = Ij ^^224 = t4, C231 = ^3) C232 = ^4; ^^244 = 

^2,^^334 = t2 +t4,C333 = 2t3, C332 = ^4, C34I = ^3^4, (^342 = ^35^444 = ^3 + 
C434 — 2i3i4. 

The ring structure (156) holds with 

W{x) = - Ux^ - hx^ + {tl - t2)x + (M4 - h) (159) 
5 
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The function 0^ is 

dW 

and the equation of the ring structure (156) holds with modW'{x) — mod[x^— 

3t3X^ - 2ts + U - t2]. 

The Landau-Ginzburg potential W in (159) is the same as the superpo- 
tential of the twisted N=2 superconformal theory for A4 type. From the 
singularity theory, this potential (159) is called a swallow tail. 

Thus we find that the random matrix theory with an external source for 
the p-th critical point gives the Landau-Ginzburg potential of the N — 2 
superconformal theory for the primary fields in the genus zero case. Our 
integral expression for the n-point correlation function may be used without 
difficulties to give the intersection numbers and the gravity descendants for 
higher genus. 

We note that these algebraic structures reduce to the Gelfand-Dikii equa- 
tion, which gives the generalized KdV hierarchies. For instance, in the case 
p=3, we obtain from our formulation the Boussinesque equation, 

and it's higher gravitational desendents. This hierarchy can be derived from 
Gelfand-Dikii equation [4]. This equation is expressed by [12] 

= (162) 

which is the generalization of Lax equation. The Q is given 

p-2 

Q = DP-Y,Mx)D' (163) 

1=0 

and the fraction power of Q is 

Q-p ^D + Y,WiD-' (164) 

i>0 

From this formulation, we obtain the relation to F as 

= -C,,^res(g"+— ) (165) 



where 



\n„n+l 



(-!)>" 

m ~ 7 TT7 T\ 7 TT (166) 

(m -|- 1) (p -|- m -|- 1) • • • [pn -\-m-\-l) 
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